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ABSTRACT 

Detailed non-LTE calculations for red supergiant stars are presented to investigate 
the influence of NLTE on the formation of atomic iron and titanium lines in the J-band. 
With their enormous brightness at J-band red supergiant stars are ideal probes of cos- 
mic abundances. Recent LTE studies have found that metallicities accurate to 0.15 dex 
can be determined from medium resolution spectroscopy of individual red supergiants 
in galaxies as distant as 10 Mpc. The non-LTE results obtained in this investigation 
support these findings. Non-LTE abundance corrections for iron are smaller than 0.05 
dex for effective temperatures between 3400K to 4200K and 0.1 dex at 4400K. For ti- 
tanium the non-LTE abundance corrections vary smoothly between -0.4 dex and -1-0.2 
dex as a function of effective temperature. For both elements, the corrections also 
depend on stellar gravity and metallicity. The physical reasons behind the non-LTE 
corrections and the consequences for extragalactic J-band abundance studies are dis- 
cussed. 

Subject headings: galaxies: abundances — line: formation — radiative transfer — stars: 
abundances — stars: late-type — supergiants 
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1. Introduction 



One of the most promising ways to constrain the theory of galaxy formation and 
evolution in a dark energy and cold dark matter dominated universe is the determination of the 
chemical composition of star forming galaxies in the nearby and high redshift universe. The 
observations of the relationship between central metallic ity and galactic mass (ILequeux et al. 



979 



Tremonti et al 



(Gamett et al, 



1997 



2004 



Maiolino et al 



Skillmai 



1998 



Gamett 



2008h a nd of metallicity gradients in spiral galaxies 



2004|) are intriguing and a clear challenge of the 



theory, which uses these observations to test the theoretical predictions of hierarchical clustering, 



galax y formation, merging, infa l 



IMF (Prantzos & Boissiei 



2000; 



2OO9I: 



, galactic winds and variability of star forrnation activity an d 



Naab & Ostriker 



2006 



Colavitti et al. 



2008 



Yin et al. 



2009; 



Sanchez-BlazQuez et al 



De Lucia et al. 



2004; 



de Rossi et al. 



2007 



Finlator & Davell2008l: 



Brooks et al 



2007 



Kdppen et al 



2007 



Wiersma et al, 



20091 : 



Dave et al. 



2011aUb ) 



So far most of our information about the metal content of star forming galaxies is obtained 
from a simplified analysis of H II region emission lines, which use the fluxes of the strongest 
forbidden lines of (most commonly) [O II] and [O III] relative to . This approach is very 
powerful and appealing, since these emission line can be observed in galaxies through the whole 
universe from low to high redshift. However, it has two weaknesses. First, this approach yields 
basically only the oxygen abundance, which is then taken, in the absence of any other information, 
as a place holder for metallicity. Secondly, th e abundances depend heavily on the calibration of the 



strong line meth od (IKewlev & Ellisor 



of up to 0.8 dex (IKudritzki et al. 



2008 



20081) and are theref ore subject to systematic uncertainties 



Bresolin et al. 



2o3). 



An altemative approach avoiding these weaknesses is the spectroscopic analysis of individual 



supergiant stars. Here, much progress has been mad e through the non-LTE low reso 



optical spectroscopy of BA blue supergiants (see e.g. 



Kudritzki et al 



2008 



Kudritzki 



ution 



2010 



Kudritzki et al. 



2OI2I . for application to spiral galaxies in the Local group). However, this 
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technique is restricted to spectroscopy in the optical. In consequence, it will not be able to take 
advantage of the fact that the next generation of extremely large telescopes such as the GMT, 
TMT, and the E-ELT will be diffraction limited telescopes at IR wavelengths allowing for adaptive 
optics (AO) supported multi-object spectroscopy. 

In this sense, the use of red supergiants (RSGs) as extragalactic abundance probes has a 
much larger potential. As direct successors of blue supergiants they have the same masses (8 to 
40 Mq at the main sequence) and the same enormous luminosities (10^ to 10^ L/L0), but their 



SEDs peak in the infrared, where the g ain in limiting 



power of telescope diameter. Recently, 



Davies et al, 



magn itude through AO goes with the fourth 



(|2010l) (hereinafter DKFIO) have introduced 



a novel technique, which uses medium resolution (R ^ 3000) spectroscopy in the J-band to derive 
stellar parameters of RSGs with the accuracy of ~ 0.15 dex per individual star. The J-band 
spectra of of RSGs are dominated by strong and isolated atomic lines of iron, titanium, silicon 
and magnesium, while the molecular lines of OH, H2O, CN, and CO which plague the H- and 
K-bands are weak and appear as a pseudo-continuum at R ~ 3000 resolution. T he prospects of the 
J-band technique in application to RSGs were discussed by Lvans et al.l (|201 IQ . who showed that 
with future instruments it would be possible to measure abundances of various chemical elements 
out to enormous distances of 70 Mpc. 

In the view of these perspectives, it is important to investigate the DKFIO method in more 
detail. One crucial question is how model atmosphere and line formation calculations are affected 
by the departures from Local Thermodynamic Equilibrium (LTE), which are expected because 
of extremely low gravities and hence low densities of RSG atmospheres. Detailed non-LTE 
line formation calculations for RSGs focussing on the key diagnostic lines in the J-band are, 
thus, a first step. In this paper, we present the first calculations of this type and investigate the 
formation of iron and titanium lines and the influence of non-LTE effects on the determination 
of abundances in the J-band. In future work, we plan to extend these studies to the formation of 
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silicon and magnesium lines, which are important for the determination of effective temperature 
in the DKF 10 method. 

The paper is structured as follows. In sections 2 and 3 we describe the model atmospheres 
and the details of the line formation calculations. Section 4 presents the results: departure 
coefficients, line profiles and equivalent widths in LTE and non-LTE and non-LTE abundance 
corrections for chemical abundance studies. Section 5 discusses the consequences for the new 
J-band diagnostic technique and aspects of future work. 



2. Model Atmospheres 

The atmospheric structure for our non-LTE line formation calculations is provided by 
the MARCS model atmosphe r es. Th e physical assumptions underlying these atmospheres are 



described in 



Gustafsson et al, 



(|2008|) . In short, these LTE models are spherically-symmetric, in 
ID hydrostatic equilibrium and include convection within the framework of the mixing-length 
theory. A carefu l discu s sion of the s trengths and the weaknesses of these models is given by 



Gustafsson et al 



that of 



(2008); 



Grevesse et al. 



(2007). 



Plea ( 20101) . The reference solar abundance mixture in these models is 



For our investigation we use a small grid of models computed assuming the mass of 15 
Mowith five effective temperatures (Teff = 3400, 3800, 4000, 4200, 4400K), three gravities 
(log^ = 1.0, 0.0, -0.5 (cgs)), three metallicities ([Z] = log Z/Z© = -0.5, 0.0, +0.5). We 
also use two values of micro-turbulence = 2 and 5 km/s, respectively. This grid covers the 
range of atmospheric parameters expected for RSGs (see DKF 10). We also compared the 1 and 
15 M0models for the same parameters finding very small differences in the NLTE abundance 
corrections. 



For the current analysis, the model atmospheres were carefully extrapolated to a continuum 
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optical depth at 500 nm logTsoo = —5, since for certain combinations of RSG atmospheric 
parameters the cores of the Ti and Fe lines have contributions from atmospheric layers around 
optical depth logTsoo ~ —3.5, which is close to the top boundary of the MARCS atmospheric 
models (logTsoo = —4.0). By this extrapolation we explored whether higher layers could 
contribute to the formation of the IR lines investigated. No significant effects were found. 

At the highest effective temperature (Teff = 4400K) and lowest gravity (log^ = —0.5) the 
MARCS models show a very slight temperature inversion, although they are converged. The 
locations of these inversions can be noticed from the kinks in the plots of departure coefficients 
(Fig. [5]and|71 see Sect. 3.1 and 3.2) for Fe land Ti I. They do not have any significant effects on 
the results. 



3. Non-LTE Line Formation 



3.1. Statistical equilibrium calculations 



For the non-LTE line formati on calculations we use D ETAIL, a non-LTE code w idely used 



for hot stars (Przvbilla et al. 



2006|) as well as for cool stars (jButler & Giddings 



1985). DETAIL 



is a well established and well tested code, which is fast a nd very efficient through the use o f the 



Rvbicki & Hummed ( 1991 



19921), and 



accelerated lambda iteration scheme in the formulation of 
allows for the self-consistent treatment of overlapping transitions and continua. It is worth noting 
that DETAIL requires the partial pressures of all atoms and important molecules to be supplied 
with a model atmosphere; these were computed using the MARCS equation-of-state package. 

Since the non-LTE effects depend crucially on the radiation field at all wavelengths, it is 
important that all relevant opacities are included in the statistical equilibrium calculations. The 
background line opacity is computed directly for each wavelength in the model atom using 
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Plezl l 



extensive linelists extracted from the NIST (|Ralchenko et al. 1120121) and K urucjjJ datab ases. For 
the current project, these tables were complemented by the TiO line data of 



19981) to account 



for the extreme line blanketing by TiO molecules. A comparison of spectral energy distributions 
calculated with DETAIL and MARCS for two typical RSG models is given in Fig. [TJ The 
agreement is satisfactory. 

The non-LTE radiative transfer in DETAIL is done in plane-parallel geometry. This might 
be a reason of concern, since the atmo spheres of RSGs have a large scale height because of 
small log^. iHeiter & Eriksson I (|2006|) have used the spherically extended MARCS models 
to compare LTE line formation calculations with radiative transfer in spherically extended 
geometry with the calculations done with the plane-parallel approximation. At most extreme cases 
abundance corrections were smaller than 0.1 dex indicating that the combination of spherically 
extended models with pl ane parallel radiative trans fer yields reasonable results. We note that the 



lowest gravities used by 



Heiter & Eriksson 



(|2006|) were log^ = 0.5 (cgs), whereas the lowest 



gravities in our grid are a factor of ten lower, i.e. logg = —0.5 (see below). On the other hand. 



Heiter & Eriksson 



(|2006l) investigated the most extreme models with very low stellar mass, i.e. 
1 Mq. However, since the spherical extend of a line forming atmosphere AR in units of the 
photo spheric radius R is: 



Mi/R^T,n/{gMfl^ (1) 

and RSGs are massive stars with masses equal and larger than 10 M©, we expect the effects of 
sphericity to be small. Indeed, compared with Heiter & Eriksson's largest extension of 12%, AR/R 
of our most extreme model (calculated with 15 Mq) is 7%. We note however that most likely 



real RSG atmospheres are more extended due to hydrodynamical phenomena (|Chiavassa et al 



^http://kurucz.harvard.edu 
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Fig. 1. — Emergent fluxes computed with DETAIL (black solid) as compared with the original 
SED's from MARCS (red, dashed). Top: Teff = 3400 K, log^ = 0, [Z] = 0. Bottom: Teff = 4400 
K, log^ = 0, [Z]=0. 



20 111 and references therein). We have also checked AR/R for the line forming region of the IR 
lines of our investigation and found values less then 2 %. This, in addition to a relatively small 
geometrical dilution factor of the radiation field ((1/1.07)^ or 14% in the most extreme model 
with AR/R = 7%), leads us to conclude that the NLTE abundance corrections are likely very 
similar in plane-parallel and spherical geometry. 



The line pro files and the NLTE abundance corrections were computed with a separate code 
SIU (lReetall999|) using the level departure coefficients from DETAIL. SIU and DETAIL share the 
same physics and line lists. 



3.2. Model atoms 



The atomic rnodels for Ti and Fe are essentially those described in iBergemannI (|201 ih and 



Bergemann et al 



(|2012|) . however some modifications were implemented for the current analysis. 
The models include the first two ionization stages of each element accounting for 332 and 397 
levels of Ti and Fe, respectively. For Fe l we have added the fine structure splitting of the a^F 
and z^D° states, since the J-band IR lines of our investigation form between these levels. We also 
removed a number of theoretically-predicted Fe I levels and radiative transitions, which according 
to our tests are not important for the statistical equilibrium of the element in RSG atmospheres. 
The Fe level diagram with transition related to the J-band lines investigated here is shown in Fig. 
[21 Similarly, the Ti I model was extended by the fine structure of all levels below 2.73 eV, and 
all transitions involving these levels were updated (Fig. |2]). The number of radiatively-allowed 



transitions is, thus, 4943 and 5328 for Ti and 



were taken from various sources (jBergemann 



e, respectively; the gf-yalues of these transitions 



2011 



Bergemann et al 



20121) . Photoionization 



cross-sections for Fe I were kindly provided by Bautista (201 1, private communication); they 
are computed on a more accurate energy mesh and provide better resolution of photoionization 
resonances compared to the older data, e.g. provided in the TOPbase. For Ti the hydrogenic 
approximation was adopted. 

To compute coUisional rates, we relied on standard recipes, commonly used in NLTE 
calculations for late-type stars, as no better alternatives are available. Electron impact excitatio n 



and io n ization cross-sec t ions w ere c omputed usin g the semi-empirical formulae from 



(119621) . 



van Regemortei] (|1962|) . and 



Seaton 



AUenl ( 1973 ). Inelastic collisions with Hi atoms were 
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computed according to the Drawin's formula in the formulation of ISteenbock & Holwegen (|1984h 
From detailed comparison with more accurate quantum-mechanical and experimental data for 
other atoms, the recipes for e~^^ colli sions are expected to give an order of magnitude accuracy 



(Mashonkina 



19961 : 



Lind et al. 



20111) . whereas the situation with Hi collisions is more uncertain. 



Recent quant um-mechanica 



exaggerated (|Barklem et al. 



calculation s suggest that for bound-bound transitions the latter are 



2011 



20121), whereas charge transfer processes can not be described 



by the Drawin's formulae at all. We note, however, that only light atoms, Li, O, Na, and Mg, have 
been investigated so far. 

Here, we performed a set of test calculations varying the collision cross-sections by several 
orders of magnitude. As expected, the most influential factor is inelastic H I collisions, which is 
also well-known from the NLTE analysis of late-type stars (Gehren et al. 2001, Bergemann et a. 
2012). However, even a factor of ten variation of these cross-sections leads to a change in the 
NLTE corrections for the IR Ti I lines by less than 0.07 dex across the RSG parameter space. A 
notable effect on the coolest models (3400 - 3800 K) is seen when the bound-bound Hi collision 
rates are decreased by four orders of magnitude. In this case, the NLTE corrections double. For 
the Fe l IR lines, the effect of varying the H I collisio ns is very small, on aver age of the order 



20121), we constrained 



0.03 dex. In our recent work on low-mass FGK stars (|Bergemann et al.l 
the efficiency of inelastic H I collisions empirically, from the analysis of stars with parameters 
determined by independent methods. We can not follow the same approach in this work, as the 
sample of nearby RSG's is very limited and stellar parameters are not accurate enough to permit 
calibration of collision efficiency using high-resolution spectra. A more accurate and realistic 
quantification of the uncertainties related to the collision rates will therefore have to await for 
detailed theoretical calculations and high-quality observations of RSG's. 







Fel 



Fig. 2. — The Fel non-LTE atomic model. The terms leading to the J-band IR transitions dis- 
cussed are highlighted in red. Only radiative transitions to and from these levels are shown in this 
plot. The calculations use a much more complex model of Fe I with many more transitions and 
simultaneously also a detailed atomic model of Fe ll which is also not shown here. 



4. Results 

The NLTE effects will be discussed using the concept of level departure coefficients bj 



= n^LTE/„LTE (2) 

where nf^^^ and are NLTE and LTE atomic level populations [cm^-^], respectively. 



Til 







5pO 5jjO 5p„ 5qO 5-^ 3gO 3pO 3-pO 3^0 3-^ IgO IpO 

¥ h° ¥ 'd -'f 'g h° 'p ' 



'D 
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'G 
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1,0 



Fig. 3. — The Til non-LTE atomic model. The terms leading to the J-band IR transitions dis- 
cussed are highlighted in red. Only radiative transitions to and from these levels are shown in this 
plot. The calculations use a much more complex model of Ti I with many more transitions and 
simultaneously also a detailed atomic model of Ti ll, which is not shown here. 

Furthermore, the importance of the NLTE effects for the determination of element abundances 
can be assessed by introducing NLTE abundance corrections Ape^xi^ where: 



AFe,Ti = logA(Fe,Ti) 



'NLTE 



logA(Fe,Ti) 



'LTE 



(3) 



is the the logarithmic correction, which has to be applied to an LTE iron or titanium 
abundance determination A of a specific line to obtain the correct value corresponding to the 
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Fig. 4. — Fine structure splitting of Fe I (left) and Ti I (right) IR J-band lines. 

use of NLTE line formation. We calculate these corrections at each point of our model grid for 
each line by matching the NLTE equivalent width through varying the Fe abundance in the LTE 
calculations. Note that from the definition of Ape^xi a NLTE abundance correction is positive, 
when for the same element abundance the NLTE line equivalent width is smaller than the LTE 
one, because it requires a lower LTE abundance to fit the NLTE equivalent width. 

It shall be kept in mind that both Ti I and Fe I have a very complex atomic structure. For 
each of the energy levels, there are about a thousand of radiative and coUisional processes, 
which contribute to the net population or de-population of a level. As a result, there is a strong 
interlocking of radiation field in different lines and continua, and it becomes highly non-trivial 
to isolate the processes explaining the populations of individual atomic levels once the statistical 
equilibrium has been established. 

Thus, in the next two subsections (4.1 and 4.2), we will give only a qualitative and rather 
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general description of the processes driving departures from LTE in the excitation-ionization 
equilibria of Ti and Fe, focusing mainly on the levels and transitions used in the spectroscopic 
J-band analysis of the RSG stars. The selected lines are given in the Table 1. Their wavelengths, 
excitation energies, and transition probabilitie s v^^ere extracted from the VALD database 



(Piskunov et al 



1995 



Rvabchikova et 



transition probabilities are taken from 



al 



1997 



Kupka et al 



O' Brian et al. 



19991, 



20001). For the Fe I lines, the 



19911) . Under typical physical conditions in 



the RSG models, all these lines are relatively strong, with equivalent widths exceeding few 
hundred mA. 



4.1. Statistical equilibrium of Fe 



Statistical equilibrium of Fe in stellar atmospheres has been extensively studied in 



application to solar-type stars, i.e., FGK spe ctral types (|Gehren et al 



2011 



Bergemann et al. 



2001 



Mashonkina et al. 



2012 Lind et all 



2012|) . However, all theses studies focused on the NLTE 
line formation of Fe I and Fe ll in the near-UV and optical spectral windows. Here we model the 
four near-IR Fe I transitions in the multiplet 296 (Table 1), which connect the lowest metastable 
2.2 eV <3^P levels to the 3.2 eV z^D° levels. The departure coefficients for these selected Fe I 
levels are shown in Fig. |5]for two selected RSG models. The LTE and NLTE unity optical depths 
in the transition at 1 1973 A (a^Pj z^^l) are indicated by the upward and downward directed 
arrows, respectively. 

In the atmospheres of the solar-type stars, the main NLTE effect on the Fe l number densities 
is over-ionization by super-thermal UV radiation field, which reduces the Fe I number densities 
compared to LTE. This effect is equally important for other minority atoms, such as Ti I which 
constitute only a tiny fraction of the total element abundance. For example, everywhere in the 
solar photosphere, the ratio of total number densities of neutral to ionized iron, Npei/Npeih is 
always below 0.1. 
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The situation is somewhat different in the cool atmospheres of RSGs. At Tgff = 3400K the 
ratio Npei/Npeii is about 10^ at optical depths logTsoo < —0.5. In consequence, photoionization 
becomes utterly unimportant for the statistical equilibrium of neutral iron at this low temperature. 
For Teff = 4400K A^Fei/A^Feii varies between 0.01 and 3.0 depending on optical depth, metallicity 
and gravity. Here, photoionization can have an effect, although the radiative rates are generally 
small due to the low T^ff of the models in combination with the extreme line blanketing (recall that 
the metallicities are close to solar). 

The analysis of the Fe I transition rates showed that for all models with Teff = 3400 to 
3800 K, irrespective of their log^ and [Fe/H], collisions fully dominate over radiatively-induced 
transitions out to the depths log T500 ~ —3 (Fig. |5l top panel). Only near the outer boundary, 
logTsoo < —3, the departure coefficients of the upper level of the multiplet 296, z^D°, deviate 
from unity, whereas the lower level a^P retains its LTE populations all over the optical depth 
scale. As a result, deviations from LTE in the opacity of the Fe I lines are negligible, and the 
NLTE and LTE line formation depths nearly coincide (as shown by the vertical marks in Fig. |5]). 
Very small NLTE effects, which show up in the line cores, are entirely due to the deviation of the 
line source function 5^ from the Planck function By{Te). Since S[j/By{Te) ~ bj/bi < 1^, that is 
driven by the photon escape in the line wings, the line cores are slightly darker under NLTE and 
the NLTE abundance corrections Ape are negative (Table 2 and 3). Metallicity and surface gravity 
have very small influence on the magnitude of NLTE effects for these cool models. 

Larger departures from LTE are found for the hottest models with Teff = 4400K (Fig. [51 
bottom panel). They are caused by radiative pumping from the lower a ^P term to higher levels 
and strong coUisional coupling. As a result a^F becomes progressively more underpopulated 
in the inner atmospheric layers and the line absorption coefficient is reduced. At the same time 

^Hereafter, i and j subscripts stand for the lower and upper level, respectively, and bi, bj are the 
corresponding departure coefficients, i.e., the ratio of NLTE to LTE occupation numbers. 
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coUisional coupling induces bi ~ bj and the line source function remains at the LTE value. As 
a consequence, the NLTE profiles become brighter in the core than the LTE profiles and the 
NLTE equivalent widths are smaller in NLTE. The magnitude of the NLTE effects decreases 
with increasing gravity and metallicity. The most extreme effects are found at Teff = 4400 K, 
logg = —0.5, and [Fe/H] = —0.5. NLTE and LTE line profiles for these parameters and for the 
model atmosphere with Tgff = 4400 K, logg = —0.5, and [Fe/H] = 0.5, are shown in Fig. |6l The 
value of the micro-turbulence has almost no influence on the size of the NLTE effects. 

4.2. Statistical equilibrium of Ti 

For Ti I, our calculations show that non-LTE effects are more significant than those for 
Fe I discussed in the previous section. The atomic structure of Ti I (ground state electronic 
configuration ls^2s^2p^3s^3p^3d^4-s^) is simpler than that of Fe I . In the NLTE model of Fe I 
, a very large number of levels allow for a more efficient colli sional and radiative coupling, 
generally leading to a stronger overall thermalization of the level populations. We also note that 
the ionization potential of Ti I is only 6.82 eV (1 . 1 eV lower than that of Fe I), thus the ionization 
balance is such that Ti ll is the dominant ionization stage for all models in our grid. The ratio of 
Nth/Ntih is 10^2 to 10"^ for the 4400K models and 0.1 to 1 for the 3400K models. 

The departure coefficients for the b^F and z^D° levels are shown in the Fig. |7] These levels 
with excitation energies of 1 .43 respectively 2.47 eV give rise to the diagnostic IR Ti I transitions 
(Table 1). The Ti I NLTE effects are a strong function of effective temperature. At Tgff = 3400K the 
lower (metastable) b^F states are overpopulated in the region of line formation through radiative 
transitions from higher levels and the upper levels z^D° are depleted by radiative transitions to 
lower levels, in particular to the ground state. In consequence, the Ti I J-band absorption lines 
become stronger in NLTE compared to LTE because of the enhanced line absorption coefficient 
and the sub-thermal source function (see Fig. |8]). 
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The situation changes towards higher effective temperatures. At Tgff = 4400K the radiation 
field has become powerful enough to deplete the lower b^F states through radiative pumping 
into higher levels. While the upper z^D° are populated from below in this way, they are also 
de-populated by radiative pumping to even higher levels, which leads to a net underpopulation. 
The strong de-population of the lower levels weakens the line absorption coefficient, while the 
source function remains close to the LTE value. As a result, the Ti I absorption lines become 
weaker in NLTE at higher temperatures (see Fig. |8}. 

4.3. Equivalent widths and non-LTE abundance corrections 

The NLTE and LTE abundance corrections and equivalent widths for the individual Ti I and 
Fe I J-band lines studied are given in Tables 2 to 5. Fig. |9] shows the NLTE abundance corrections 
for some of these lines computed using ^ = 2 km s^^ as a function of effective temperature, 
gravity and metallicity. The results for microturbulence 5 km s ^ are not shown because they 
are nearly identical. As a consequence of the small Fe I NLTE effects discussed in section 4.1 
the NLTE abundance corrections are very small and reach a maximum value ~ 0.1 dex only 
at the highest effective temperature. The medium resolution J-band fitting method of RSG 
spectra developed by DKFIO achieves an accuracy of 0.15 dex on average for the metallicity of 
an individual RSG. Thus, we do not expect NLTE effects in Fe I to heavily affect the results. 

On the other hand, the NLTE effects for Ti I are more pronounced as discussed in section 3.2. 
This is clearly reflected in the NLTE abundance corrections, which become ~ —0.4 dex at low 
effective temperatures and change to ~ 0.2 dex at high effective temperature. Thus, if the Ti I 
lines were the only features to derive stellar metallicities from J-band spectra, NLTE effects would 
imply large corrections. However, the technique introduced by DKFIO uses the full information of 
many J-band lines from different atomic species (7 Fe I, 2 Mg I, 2 Ti I, 3 Si I lines) simultaneously 
to determine metallicity together with stellar parameters. Since the Fe I lines are only weakly 
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affected by NLTE, we expect Ti I NLTE effects to have a smaller influence on the determination 
of the overall metallicity than the NLTE abundance corrections for the Ti I near-IR lines only. 



4.4. J-band medium resolution spectral analysis 

In order to assess the influence of the NLTE effects on the J-band medium resolution 
metallicity studies, we have carried out the following experiment. We calculated complete 
synthetic J-band spectra with MARCS model atmospheres and LTE opacities for all spectral lines 
except the Fe I and Ti l for which we used our NLTE calculations. We then used these synthetic 
spectra calculated for Teff from 3400K to 4400K with different log g and different metallicity (and 
with added Gaussian noise corresponding to S/N of 200) as input for the DKFIO analysis using 
MARCS model spectra calculated completely in LTE. The metallicity grid for the MARCS model 
spectra had a resolution of 0.25 dex b etween [Z] = -1 an d -1-0.5 and a cubic spline interpolation 



was used between the grid points (see 



Evans et al. 



2011 ). From the metallicities recovered in this 



way we can estimate the possible systematic errors when relying on a complete LTE fit of RSG 
J-band spectra. 

We found small average metallicity corrections of —0.15 dex at low temperatures and +0.1 
dex at high temperatures with only a weak dependence on input metallicity (Fig. 11). This 
qualitative behaviour is as expected given the results shown in Fig. 9 and 10. Quantitatively, the 
effects are small, since the results of these tests are dominated by the Fe l lines, which are more 
numerous than the Ti I lines and for which the behavior NLTE corrections are minor (note that 
in this paper we have only discussed the non-LTE effects of the four strongest Fe I lines). The 
[Z] resolution of the model grid for this experiment is 0.25 dex and we know that the systematics 
uncertainties of the analysis technique are of the order of 0.15 dex (DKFIO), thus, the NLTE 
corrections are marginally within the accuracy of the fitting procedure. However, very obviously 
the inclusion of NLTE effects in the calculation of J-band Fe I and, in particular, Ti I will improve 
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the accuracy of future work. 



4.5. Modeling uncertainties 



Spectroscopic 



(IGustafsson & Plez 



analys is of red supergiants is generally a highly non-trivial matter 



1992|) . as other complexities arise. Among the most severe is sub- 



photospheric convection. At present, it is not possible to systematically evaluate the effect of 
ID hydrostatic equilibrium approximation on the determination of stellar parameters from RSG 
spectra. The only realistic non-gray 3 D radiative hydrodyna mics simulation of convection for a 



single RSG model was performed by 



Chiavassa et al. 



(|201l|) . They demonstrated that there is 



no unique ID static model, which would provide similar observable characteristics, i.e. SED 
shape and colors, as the analogous 3D RHD model. To match the optical spectrum dominated 
by TiO opacity, a ~ 100 K cooler ID model is needed, whereas the IR spectral window can 
be best matched using a ~ 200 K hotter ID model. The NLTE effects, however, are sensitive 
to the radiation field across the whole spectrum, from UV to the IR. Thus, it is not possible to 
predict the changes in the NLTE abundance corrections simply by using a cooler or hotter ID 
static model. One could, in princip le, expect that, since the mean T(t) relation of a typical 3D 



RHD model (Chiavassa et al 



2011 



their Fig. 5) is shallower than that of a corresponding ID 



static model, the NLTE effects would be smaller. Our recent work on low-mass solar-like stars 



(IBergemann et al 



2012|) showed that the largest differences in terms of NLTE effects occur at 



low metallicities, where 3D models are up to 1000 K cooler than ID hydrostatic models. The 
RSG models investigated in this work have metallicities close to solar. In any case, full NLTE 
calculations adopting, at least, a mean structure of the 3D RHD model, are necessary. Given our 
recent experience with calculations of such type for the solar- like stars (Bergemann et al. 2012), 
we are now in position to perform similar analysis for RSG's. 
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5. Conclusions and Future Work 

We have presented detailed calculations of the non-LTE effects in Fe I and Ti l lines in the 
J-band, relevant to the analysis technique of DFKIO. The non-LTE corrections to abundances 
measured from individual Fe I lines are small - less than 0.15 dex across the parameter range 
studied. On the other hand, for Ti I these corrections are larger, 0.4 dex at the extreme. However, 
as there are only two Ti I lines in the J-band spectral window compared to many Fe I lines and 
other lines of other elements such as Mg I and Si I the impact on the derived metallicity of a given 
spectrum is small. Differential analysis of non-LTE corrected and uncorrected spectra shows that 
the average difference in measured metallicity [Z] is < 0.15 dex, i.e. still within the accuracy range 
of the medium resolution analysis procedure but large enough and showing systematic trends with 
Tgff and [Z] that it seems worthwhile to include the non-LTE effects in the calculation of the J-band 
synthetic spectra. This will ultimately improve the accuracy of the method, in particular, when it 
will be applied not only to measure overall metallicity but also individual element abundances and 
the diagnostically important ratio of a to iron elements. 

In terms of future work, we will next turn our attention to the other prominent species in 
the y-band spectral window, namely those of Si I and Mg I. These lines also contain important 
metallicity information, since they help constrain the measurement of a star's a-element 
abundance. In order to determine the non-LTE corrections to the 7-band transitions of these 
species we have already begun to construct model atoms, and the computation of the NLTE line 
strengths will begin shortly. 

This work was supported by the National Science Foundation under grant AST-1 108906 
to RPK. Moreover, RPK acknowledges support by the Alexander- von-Humboldt Foundation 
and the hospitality of the Max-Planck-Institute for Astrophysics in Garching and the University 
Observatory Munich, where part of this work was carried out. 
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Fig. 5. — The NLTE departure coefficients of Fe I for RSG models with log g = -0.5, [Z] = -0.5 
and Tgff = 3400K (top) and 4400K (bottom) as a function of optical depth. Blue: Fe II ground 
state a^Dg^2- Black: lower fine structure levels of J-band IR transitions, a^Pj (diamonds, dashed), 
a^P2 (dotted), a ^Pg (solid). Red: upper fine structure levels of IR- transitions, z^Dq (crosses), z^Dg 
(squares, dashed), z^D4 (triangles). The LTE and NLTE line core optical depths log t(1 1973 A, 
Fe i) = are also indicated by the upward and downward directed arrows, respectively. 
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1.19715 1.19730 1.19745 

Fig. 6. — The NLTE (black, solid) and LTE (red, dashed) profiles of the Fe I line at 11973 A 
computed for Teff = 4400K, log g = -0.5 and [Z] = -0.5 with microturbulence £,t — 5 km/s. 
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log X (500 nm) 

Fig. 7. — The NLTE departure coefficients of Ti I for RSG models with log g = -0.5, [Z] = -0.5 
and Teff = 3400K (top) and 4400K (bottom). Blue: Till ground state a'^F. Black: lower fine 
structure levels of J-band IR transitions, b^F^ (dotted), b^F2 (solid). Red: upper fine structure 
levels of IR- transitions, z^Dj (triangles), z^Dj (diamonds). The LTE and NLTE line center optical 
depths log t( 1 1949.58 A, Ti l) = are also indicated by the upward and downward directed arrows, 
respectively. 
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Fig. 8.— The NLTE (black, solid) and LTE (red, dashed) profiles of the Tilline at 11949.58 A 
computed for log g = -0.5 and [Z] = -0.5 with microturbulence i^t = 5 km/s and Tgff = 3400K (top) 
and 4400K (bottom). 




Fig. 9. — NLTE abundance corrections as a function of effective temperature for microturbulence 
= 2 km/s for Fe il 1593 A(top),Fe il 1638 A(2nd row), Ti il 1893 A(3rd row) and Ti il 1949 (bot- 
tom). Left column: log g = -0.5, middle column: log g = 0.0, right column: log g = 1.0. Black 
solid: [Z] = -0.5, blue dotted: [Z] = 0.0, red dashed: [Z] = -1-0.5. 
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Fig. 10. — Influence of the Til and Fel non-LTE effects on the DKFIO J-band metallicity 
determination as a function of effective temperature. The numerical experiment is described in the 
text. Left: log g = —0.5, middle: log g = 0.0, right: log g = 1.0. Circles: [Z] = —0.5, triangles: 
[Z] = 0.0, squares: [Z] = 0.5 
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Table 1 . J-band Fe I and Ti I lines 



Elem. 


X 




lower 




upper 


log gf 




A 


[eV] 


level 


[eV] 


level 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 



Til 


11892.85 


1.43 




2.47 




-1.908 




11949.58 


1.44 




2.48 




-1.760 


Eel 


11638.25 


2.18 




3.25 




-2.214 




11973.01 


2.18 


a'P, 


3.22 


z^Dl 


-1.483 




11882.80 


2.20 


a'Pj 


3.24 




-1.668 




11593.55 


2.22 




3.29 




-2.448 
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Table 2. Non-LTE abundance corrections for the Ti I and Fe I lines =2 kms ^) 



Tcff log g [Z] Atii Ath AfcI Apel AfcI Apel 

11892.85 11949.58 11593.55 11638.25 11882.80 11973.01 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 



4400. 


-0.50 


0.50 


0.26 


0.26 


0.02 


0.02 


0.02 


0.02 


4400. 


-0.50 


0.00 


0.22 


0.22 


0.05 


0.04 


0.03 


0.03 


4400. 


-0.50 


-0.50 


0.15 


0.16 


0.12 


0.10 


0.08 


0.07 


4400. 


0.00 


0.50 


0.19 


0.18 


0.01 


0.01 


0.01 


0.01 


4400. 


0.00 


0.00 


0.17 


0.16 


0.05 


0.03 


0.02 


0.02 


4400. 


0.00 


-0.50 


0.12 


0.12 


0.11 


0.09 


0.07 


0.06 


4400. 


1.00 


0.50 


0.08 


0.06 


0.01 


0.00 


-0.01 


-0.01 


4400. 


1.00 


0.00 


0.09 


0.08 


0.03 


0.02 


0.01 


0.01 


4400. 


1.00 


-0.50 


0.08 


0.07 


0.08 


0.08 


0.05 


0.05 



Note. — This table is published in its entirety in the electronic edition of ApJ. A portion is shown here 
for guidance regarding its form and content. 



3. Non-LTE abundance corrections for the Ti I and Fe I lines = 5 kms 



Tcff log g [Z] Atii Ath AfcI Apel AfcI Apel 

11892.85 11949.58 11593.55 11638.25 11882.80 11973.01 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 



4400. 


-0.50 


0.00 


0.25 


0.25 


0.03 


0.03 


0.03 


0.03 


4400. 


-0.50 


0.00 


0.22 


0.22 


0.06 


0.06 


0.04 


0.04 


4400. 


-0.50 


-0.50 


0.17 


0.16 


0.12 


0.12 


0.11 


0.11 


4400. 


0.00 


0.50 


0.20 


0.20 


0.02 


0.01 


0.00 


0.01 


4400. 


0.00 


0.00 


0.17 


0.17 


0.06 


0.05 


0.03 


0.03 


4400. 


0.00 


-0.50 


0.15 


0.14 


0.11 


0.11 


0.10 


0.10 


4400. 


1.00 


0.50 


0.11 


0.10 


0.01 


0.00 


-0.02 


-0.02 


4400. 


1.00 


0.00 


0.11 


0.11 


0.04 


0.03 


0.02 


0.01 


4400. 


1.00 


-0.50 


0.12 


0.10 


0.08 


0.09 


0.08 


0.08 



Note. — This table is published in its entirety in the electronic edition of ApJ. A portion is shown here 
for guidance regarding its form and content. 
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Table 4. Equivalent widths ^ of the Ti I and Fe I lines =2 kms ^) 



Teff 




[Z] 


Wa,t,i 


Wa.t,i 


Wa.t,i 


Wa.t,! 


WA.FeI 


WA.Fel 


WA.Fd 


WA.Fel 


WA.Fel 


WA.FeI 


WA.Fd 


WA.FeI 








11892 


11892 


11949 


11949 


11593 


11593 


11638 


11638 


11882 


11882 


11882 


11973 








LTE 


NLTE 


LTE 


NLTE 


LTE 


NLTE 


LTE 


NLTE 


LTE 


NLTE 


LTE 


NLTE 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


4400. 


-0.50 


0.50 


221.6 


186.3 


239.7 


205.6 


438.5 


436.0 


492.1 


488.2 


658.1 


649.8 


773.3 


759.5 


4400. 


-0.50 


0.00 


155.8 


123.2 


176.2 


143.3 


386.5 


380.7 


425.9 


420.4 


530.0 


522.5 


597.4 


586.5 


4400. 


-0.50 


-0.50 


80.1 


61.4 


98.6 


77.8 


336.4 


324.9 


370.2 


358.9 


443.0 


431.8 


483.4 


470.6 


4400. 


0.00 


0.50 


218.6 


194.5 


236.1 


213.6 


428.4 


426.6 


482.0 


480.2 


647.6 


643.7 


762.2 


754.5 


4400. 


0.00 


0.00 


149.4 


126.2 


168.7 


146.1 


373.7 


368.4 


413.3 


408.8 


518.2 


512.6 


585.8 


577.6 



''equivalent widths Wa are given in mA 

Note. — This table is published in its entirety in the electronic edition of ApJ. A portion is shown here for guidance regarding its form and content. 
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Table 5. Equivalent widths ^ of the Ti I and Fe I lines = 5 kms ^) 



Teff 




[Z] 


Wa,t,i 


Wa.t,i 


Wa.t,i 


Wa.t,! 


WA.FeI 


WA.Fel 


WA.Fd 


WA.Fel 


Wa.foI 




Wa.foI 


Wam 








11892 


11892 


11949 


11949 


11593 


11593 


11638 


11638 


11882 


11882 


11882 


11973 








LTE 


NLTE 


LTE 


NLTE 


LTE 


NLTE 


LTE 


NLTE 


LTE 


NLTE 


LTE 


NLTE 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


4400. 


-0.50 


0.50 


387.6 


304.7 


432.6 


349.8 


871.8 


866.0 


944.3 


938.2 


1111.8 


1106.5 


1208.6 


1198.9 


4400. 


-0.50 


0.00 


234.0 


170.0 


278.6 


209.9 


781.9 


768.3 


850.8 


838.4 


990.4 


980.0 


1060.6 


1047.5 


4400. 


-0.50 


-0.50 


96.0 


68.2 


123.1 


91.3 


678.3 


651.4 


748.6 


722.2 


877.9 


853.4 


936.3 


910.1 


4400. 


0.00 


0.50 


385.6 


324.5 


428.4 


369.3 


847.7 


843.4 


920.2 


917.0 


1088.4 


1087.5 


1186.2 


1181.5 


4400. 


0.00 


0.00 


225.4 


176.5 


268.0 


216.7 


751.8 


739.6 


820.4 


809.8 


960.5 


952.4 


1031.2 


1020.8 



''equivalent widths Wa are given in mA 

Note. — This table is published in its entirety in the electronic edition of ApJ. A portion is shown here for guidance regarding its form and content. 



